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Abstract. The propagation of Galactic Cosmic Ray nuclei having energies between 
100 MeV/nuc and several PeV/nuc is strongly believed to be of diffusive nature. 
The particles emitted by a source located in the disk do not pervade the whole 
Galaxy, but are rather confined to a smaller region whose spatial extension is 
Oh' related to the height of the diffusive halo, t he Galactic wind and the spallation 

rate. Following the pioneering work of Jones ( 1978), this paper presents a general 
| study on the spatial origin of cosmic rays, with a particular attention to the role of 

spallations and Galactic wind. This question is different, and to a certain extent 



X 



disconnected, from that of the origin of cosmic rays. We find the regions of the 
disk from which a given fraction of cosmic rays detected in the Solar neighborhood 
were emitted (/-surfaces). After a general study, we apply the results to a realistic 
source distribution, with the propagation parameters obtained in our previous 



system atic analysis of the observed secondary-to-primary ratios IIMaurin et al 



2002a). The shape and size of these /-surfaces depend on the species as well as 
on the values of the propagation parameters. For some of the models preferred by 
our previous analysis (i.e. large diffusion slope 8), these /-surfaces are small and 
in some extreme cases only a fraction of a percent of the whole Galactic sources 
actually contribute to the Solar neighborhood Cosmic Ray flux. Moreover, a very 
small number of sources may be responsible for more than 15 % of the flux detected 
in the Solar neighborhood. This may point towards the necessity to go beyond 
the approximations of both homogeneity and stationarity. Finally, the observed 
primary composition is dominated by sources within a few kpc. 



Key words. Cosmic rays 



2 Taillet & Maurin: Spatial Origin of Cosmic Rays in Diffusion Models 

1. Introduction 



The propagation of charged Cosmic Ray nuclei, in the energy range going from a few 
100 MeV/nuc and a few PeV/nuc, is strongly affected by the Galactic magnetic field. It 
is a diffusive process, so that the Cosmic Rays emitted by a single source spread out in 
time, pervade the whole Galaxy, and can escape the Galaxy when reaching its boundaries. 
Those coming from a source located far from the Sun have a larger probability of escaping 
than reaching the Solar neighborhood. It is the opposite for nearby sources, so that 
the cosmic ray fluxes in the solar neighborhood are more sensitive to the properties of 
the local sources (as opposed to the remote sources). Other effects like spallations and 
Galactic wind further limit the distance Cosmic R ays travel bef ore being detected. Some 



■Tone; 



19781) where convective escape 



consequences of the Galactic wind were studied in 
was compared to escape through the top and bottom boundaries of the Galaxy. 

The goal of this paper is to go one step beyond by providing a general study on 
the spatial origin of cosmic rays, i.e. to answer the question "from which region of the 
Galaxy were emitted the cosmic rays detected in the solar neighborhood?" . This question 
is different, and to a certain extent disconnected, from that of the origin of cosmic rays 
("What are the astrophysical objects which are responsible for the acceleration of cosmic 
rays?") which is still much debated. We believe that it is nevertheless an interesting 
question, for several reasons. First, we find that the answer may cast some doubt on the 
validity of the stationary model, upon which most studies on Cosmic Rays are based. 
Second, it gives some clues about the spatial range beyond which the cosmic ray studies 
are blind to the sources. Finally, this study may be of interest to optimize the propagation 
codes based on Monte-Carlo methods, by focusing the numerical effort on the sources 
that really contribute to the detected flux. 

The reader who does not want to go through the pedagogical progression can go 
directly from the general presentation of the method in §[21 to its application in realistic 
cases in § For the others, the effect of escape is studied in § and that of spallations 
and Galactic wind is studied in § \5\ Then, § \G\ studies the effect of a realistic source 
distribution. Finally, the fully realistic case is considered in § El The results and the 
perspectives are discussed in the last section. For convenience, we will use the word 
/-surfaces to describe the surfaces in the thin disk within which the sources form the 
fraction / of cosmic rays detected at the observer location. 



2. Description of the method 

A stationary point source emits particles that diffuse in a given volume. At the boundaries 
of this volume, the particles are free to escape and the density drops to zero. After 
a sufficiently long time, the stationary regime is eventually reached and the density 
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profile is established inside the diffusive volume. If several sources are present (or even a 
continuous distribution of sources), their contributions add linearly at each point. 

The question we wish to answer is the following: a Cosmic Ray being detected at the 
position r a of an observer (in practice, this will be the position of the Sun, and we refer 
to this position as the Solar neighborhood), what is the probability density 

dP {emitted : r s , r s +dr s | observed: r } _ dV {r s \r a } 
dr s ~ dr s 

that this Cosmic Ray was emitted from a source located at the position r s ? Such a ques- 
tion falls among classical problems of statistics. A rigorous theoretical frame is provided 
by the Bayes approach that summarizes the proper use of conditional probabilities. A 
cruder but sufficient (and equivalent) treatment is given by the frequency interpretation. 
The probability written above is simply given by 

dV{r s \r } dN[r s ^r ] /dr s 

dr s N [-> r a ] ' U 

where N [— > r ] is the number of paths reaching r and dM [r 3 — > r ] /dr s is the density 

of paths going from r s to r . We finally notice that the latter number determines the 

density of Cosmic Rays that reach the position r D , when a source is placed at r s . We 

can thus write 

where the density N Ys (r Q ) is the solution of the propagation equation for a point source 
located at r s . The normalization factor in this relation is obtained by imposing that 
dV/dr s actually is a probability density, i.e. is normalized to unity. We refer to the 
contours on which the probability density is constant as isodensity contours. 

If the sources are distributed according to w(r s ), the probability that a Cosmic Ray 
detected at r Q was emitted from a surface S is given by 

This probability contains all the physical information about the spatial origin of Cosmic 
Rays. We define the /-surfaces, inside which the sources contribute to the fraction / of the 
detected flux, by the relation V {<S|r } = /. Actually, even for a given value of /, there are 
many different surfaces, delimited by different closed contours, fulfilling this condition. 
We focus on the smallest of these surfaces, which is precisely delimited by an isodensity 
contour. We also use the term rii m -probability for the quantity V {r s < ru m \r }. 

3. The escape through the diffusive volume boundaries 

The region in which diffusion occurs is limited by surfaces (hereafter the boundaries) 
beyond which diffusion becomes inefficient at trapping the particles, so that they can 
freely escape at a velocity close to c. The density outside the diffusive volume is very 
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small, and it is very reasonable to suppose that the boundaries are absorbers, i.e. they 
impose a null density (TV = 0). 

It is well-known that the shape and location of the boundaries play a crucial role for 
diffusive propagation. This section shows that the Cosmic Rays emitted from standard 
sources in the disk are not sensitive to the radial extension of the Galaxy, but only to its 
top and bottom edge. To this aim, it is sufficient to concentrate on pure diffusion and 
to neglect spallations, the Galactic wind and reacceleration. Indeed this is a conservative 
case as these effects can only make the diffusion process even less sensitive to the presence 
of the boundaries (see below) . Moreover, we consider the case of a homogeneous source 
distribution located in the disk w(r s ) oc S(z), which also leads to a conservative result if 
compared to a realistic radial distribution of sources. 

We first consider the pure diffusion equation with a Dirac source term 



In unbounded space, the solution is given by N rg (r a ) = l/4TrK\\r — r s \\. The influence 
of the boundaries is estimated by solving this equation in three situations: first we con- 
sider only a side boundary, then only a top plus bottom boundary, and finally all the 
boundaries. 

3.1. Boundaries influence 



Fig. 1. Geometry of the diffusive volume. 

Our Galaxy can be represented as a cylindrical box with radial extension R and 
height L (see App. A for further details). The probability density dV cy i(r s \r )/dr s can 
be computed for arbitrary source and observer positions r s and r Q , using a Fourier- 
Bessel decomposition of the density. In our case, the observer is located near the Sun, at 
a Galactocentric distance i? Q ~ 8.5 kpc. Unless the diffusive halo height is very large, the 
top and bottom boundaries located at z = ±L are nearer to us than the side boundary 
located at R — 20 kpc. As a result, we expect the effect of the side boundary to be 
smaller. The first simplified situation we consider is that of an observer located at the 



- KAN(r) = S(r - r s ) . 



(5) 




L 



Taillet & Maurin: Spatial Origin of Cosmic Rays in Diffusion Models 5 

center of the Galaxy. (In the case of an infinite disk, i.e. R — > oo, this amounts to a mere 
redefinition of the origin of the disk) . 

With r D = 0, the solution for a point source in this particular geometry is given in 
App. A. The probability density that a particle reaching the observer was emitted from 
a point located at a distance r 8 from the center is thus given by (with p s = r s /R) 

normalization being obtained by imposing J Q tfP(r s |0) = 1. The ri; m -probability is given 

by 

t> < ^ \n\ riim Jv^ J i(Cinim/-R) ^ . ,,. r/m ) f^tanh(CiL/-R)\ 
^ < n im |0) = — jg c2ji2(G) x tanh^L/i?) j . |E j -(7) 

This probability is independent of the value of the diffusion coefficient K. 
3.1.1. Side boundary 

I 
I 

i i i 

i i i 

i . i 

i i 




i 
i 



Fig. 2. Geometry of the diffusive volume in the limit L — > oo. 

The escape from the side boundary (located at r = R) is disentangled from the 
escape from the z — ±L boundary by first considering the limit L —t oo. For the sake 
of simplicity, we will, as above, only study the effect of this boundary on observations 
performed at the center of the Galaxy. In the limit L — ► oo, we have coth(QL/ R) f» 1 in 
expression lfH|. This gives for the Him-probability, 

-77 I ^ irrt riim -M^W-g) ] fy> 1 1 

^ < r nm \0) = — (E ] • {g cfjUG) } • 
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Fig. 3. Geometry of the diffusive volume in the limit R — > oo. 
3.1.2. Top and bottom boundaries 

The influence of the z — ±£ boundaries, in the case of an infinite disk (R — ► oo) is now 
considered. In this limit, the sum over Bessel functions can be replaced by an integral 
and we obtain (see App. B.3) 



dV L (r s \0) 



cPr, 



Jul./') tanh ( ) (/./■ 



(8) 



which allows to compute the ni m -probability V\^{r s < nim|0) as before, which is a func- 
tion of rii m /L only. These integrals are somewhat intricate to compute numerically, due 
to the very slow convergence. In this particular case, the accuracy of the numerical cal- 
culation can be checked for r > L, as a detailed study of the function iJHJ shows that in 
this limit 

f 



x 2, 



-nr s /2L 



4irKr s \ L 
It is also noticeable that the quantity 
N L 



/esc(^*s) — 1 



J J (x) jl - tanh (~^-) j dx 



(9) 



(10) 



gives the fraction of Cosmic Rays emitted from a distance r s that has escaped the diffusive 
halo before reaching us. 



3.2. Summary: the effect of boundaries on primary species 

Fig. CD shows the probability density computed above as a function of r s for unbounded 
space, for the cylindrical geometry with several halo sizes L, i.e. Eq. J7J), and for the 
two limiting cases corresponding to L — > oo or R — > oo. We also show, in Tab. ^ the 
radii of the /-surfaces, in the two cases R = 20 kpc and R — > oo. It can be noticed that 
even if the source distribution is infinite in extent, the finite size of the halo limits the 
quantity of Cosmic Rays that reach a given point. The mean distance from which the 
Cosmic Rays reaching the center is given by (r s ) — 1.4 L. This effect dominates over 
the leakage through the side boundary r s = R, and it will be even more negligible in 
realistic situations, as (i) the source density is small near the edge of the disk, (ii) when 
the spallations and Galactic wind are considered, most Cosmic Rays are destroyed or 
blown out of the disk before they have a chance to reach this side boundary. 
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* x * x dPoo(r s )= 1 /(4Tir s 

dP S p h (r s ) = (1-p s )/(47Tr s 




Fig. 4. Cosmic ray probability density as a function of r s (distance of the S(r 3 ) source 
in the disk), for several values of L and for a disk of radius R = 20 kpc. Big stars are for 
unbounded model, dotted line is for a spherical boundary at radius R, small stars are for 
top and bottom boundaries, and solid lines are for cylindrical boundaries. 





/(mm) = 50% 
R = oo R = 20 kpc 


f(r Um ) = 90% 
R = oo R = 20 kpc 


/(r, im ) = 99% 
R = oo R = 20 kpc 


L = oo 
L = 20 
L = 5 kpc 
L = 1 kpc 


6.2 kpc 
12.6 kpc 6.1 kpc 
3.1 kpc 2.95 kpc 
0.63 kpc 0.63 kpc 


14.1 kpc 
39 kpc 14 kpc 
9.5 kpc 8.6 kpc 
1.9 kpc 1.9 kpc 


18.2 kpc 
68 kpc 18.2 kpc 
17 kpc 14.6 kpc 
3.4 kpc 3.4 kpc 



Table 1. This table indicates the radius ru m inside which a given fraction /(n™) = 
V{r s < riim | O) of Cosmic Rays reaching the center were emitted from, for several L and 
in the case of the infinite disk and R = 20 kpc. 



An important consequence is that as long as the observer and the sources are not 
too close to the side boundary, the density only depends on the relative distance to 
the source in the disk, so that it may be assumed, for numerical convenience, that the 
observer is either at the center of a finite disk, or in an infinite disk. In all the paper, 
i.e. for standard sources in the disk, we will consider the limit R — > oo, i.e. we use the 
integral representation described in App. B.3. 
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4. Secondary and radioactive species 

4.1. Progenitors of stable secondaries 

As can be seen in App. A. 3, the secondary distribution from point-like primary sources 
is related very simply to the primary distribution itself. One could find strange to speak 
about secondaries as we have not, for the moment, included spallations in the model. The 
right picture is the following: a primary emitted at r s propagates and from time to time 
crosses the disk (mostly filled with hydrogen, density Uism)- During this crossing, there 
is a probability ^iSM-w.cr P rim->sec to create a secondary, that in turn propagates in the 
diffusive volume until it reaches (or not) the experimental setup. This will be taken into 
account properly in the next section. However, in order to have a compact expression, 
a crude estimation can be obtained by neglecting the influence of spallations on the 
primary and secondary component. This is obtained if one discards ri nc i in the terms 
A? rim and ^4| ec of Eq. i|A.5|) . The net result will be an overestimation of the distance the 
secondaries come from since their destruction is discarded two times; once under their 
primeval primary form and once in their secondary form. 

We find, in the case R — > oo (see App. B.3), and for a homogeneous distribution of 
sources, 

dV (r s \R e ) g r^M xtanha (*L\ dx 
d 2 r s Jo x V r s J 

The resulting integrated probabilities are shown in Tab. |3 The source of the primary 
that will give the secondaries observed at a given point is located farther away than the 
sources of the primary we detect (compare Tabs. [5] and This may be of importance 
if for instance the source composition or the source intensity varies with position: in 
the ubiquitous secondary-to-primary ratio, the numerator is sensitive to sources located 
on a greater range than the denominator. Moreover, these secondaries set the size of 
an effective "local" zone outside of which the particles reaching the Solar neighborhood 
have never been. The local observations tell nothing about the propagation conditions 
outside of this zone. One could object that this conclusion is mainly based on the /- 
surfaces which refer to the sources contributing to observed CR, but that the Cosmic 
Rays reaching us from these sources actually sample (via random walk) a much larger 
volume. This is actually not the case, as a particle wandering too far has a very small 
probability to ever come back to us. This point can be made more quantitative, as a 
simple reasoning shows that the probability V\ACB\ that a particle emitted in A and 
reaching B has passed through C is given by V[AC]V[CB], which is closely related to 
Na^c) x Nb{tc)- This later quantity is small as soon as C is too far from A or B. 
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/(Him) = 50% 

i? = oo i? = 20 kpc 


/(rum) = 90% 
R = oo i? = 20 kpc 


/(riim) = 99% 
R = oo i? = 20 kpc 


L = oo 
L = 5 kpc 
L = 1 kpc 


8.6 kpc 
5.5 kpc 5.3 kpc 
1.1 kpc 1.1 kpc 


15.3 kpc 
12.5 kpc 12 kpc 
2.5 kpc 2.5 kpc 


18.5 kpc 
25 kpc 17.2 kpc 
4.4 kpc 4.4 kpc 



Table 2. This table indicates the radius inside which a given fraction /(rn m ) of secondary 
Cosmic Rays reaching the center were emitted from, for several L and in the case of a 
disk of radius R — 20 kpc. The last line shows that for small L, the effect of the side 
boundary is completely negligible. 



4.2. Radioactive secondaries 

In the case of an unstable species with a lifetime t, formula (|A.4(I can be written as 
dV taA {r s \0} oc dp s Y MM (n) 



where r rfl r i 



( Lebedev 



1972) 



t 1 = 7 Vq . This expression can be transformed using the identity 



1 f°° e~ ap f 1 e~ ap 

-pJo(CiP) d P~ / pJo{(iP)dp 



y/Q + a 2 Jo P Jo P 

The approximation in the last step is valid if the exponential term decreases with p 
fast enough (i.e. a is large so that the upper limit can be set to 1 in the integral). We 
then recognize in (|llfl the Fourier-Bessel transform of exp(— ap)/p, so that finally the 
normalized probability reads 

dV lad {r s \0} = ex P(-^A-ad) d 2 Ts (12) 
Zir r s . « rad 

where the following typical length has been introduced 



l lad = \ l —— = 0.17 kpc x J T J J" . (13) 

W r rad y V 0.03 kpc 2 Myi- 1 V 1 Myr V 1 

Indeed, this result can be derived much more straightforwardly starting from the sta- 
tionary equation —KA r N(r) + r ra d^V(T') = (with a source at the origin) in unbounded 
s pace. This is also in fu ll agreement with the expression given in App. B (see also Sec. 4.1) 
of Don ato et alJ l|2002f) . where we found the same expression starting from the propagator 
of the non-stationary diffusion equation in unbounded space. 

To sum up, Eq. (|12fl is valid as long as Z ra( j <C R and / ra( j <C L: the propagation of 
the unstable species can be then considered as local, with a typical scale ^ ra d- This is no 
longer the case if the lifetime r = "/tq is large, which is the case at high energy because 
of the relativistic factor 7, even if the proper lifetime to is short. The rii m -probability is 
straightforwardly derived. As on these typical scales, the source distribution can safely 
taken to be constant, the distance rn m is expressed as 

rum = -kad x ln(l - /) . (14) 
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It means that the sources that contribute to the fraction / = (50 — 90 — 99)% of the 
radioactive species measured flux are located inside the disk of radius n; m = (0.7 — 2.3 — 
4.6) x Z rac i. The effect of a local underdensity around the Sun is discussed later. 



4.3. Electrons and positrons 

Cosmic ray sources also emit electrons and positrons. In contrast with the nuclei, these 
particles are light, so that they are subject to much stronger energy losses, due to syn- 
chro tron radiation an d inver se Compton. This results in an effective lifetime given by 
(e-g- 



Aharonian et al 



(1995|) n 



300 Myr x (1 GeV/E). The results given in the 



previous section on radioactive species can be applied to this case, with a scale length 



^loss 



1 kpc x 



1 GeV 



K 



E 



0.03 kpc 2 Myr" 1 



Formu lae 11211 and (I14|) can be used with l rai ± <-> r\ oss . This effect is discussed by 



Aharonian et al 



1 1995J) to show that a nearby source may be necessary to explain the 



high energy electron flux observed on the Solar neighborhood. 



4.4. Summary: pure diffusive regime, an upper limit 

The important conclusions at this point are that i) most of the stable primary Cosmic 
Rays that reach the Solar neighborhood were emitted from disk sources located within 
a distance of the order of L, such that the R = 20 kpc boundary can reasonably be 
discarded ii) the secondary species composition is determined by sources located farther 
away than those determining the primary composition; iii) radioactive species may come 
from very close if their lifetime is so short that \JK^tq < L, high energy electrons and 
positrons definitely do. 

These conclusions are expected to be stronger when spallations, Galactic wind and a 
realistic source distribution arc taken into account. All these effects will limit even more 
the range that the particles can travel before reaching the Solar neighborhood. 

5. The effects of spallation and convection 

5.1. Pure convection 

The diffusion of Cosmic Rays may be disturbed by the presence of a convective wind of 
magnitude V C1 directed outwards from the disk. For numerical convenience, a constant 
wind has been considered, although other possibilities (especially a linear dependence) are 
probably more justified on theoretical grounds (see discussion in MTD02). The effect is to 
blow the particles away from the disk, so that those detected in the Solar neighborhood 
come from closer sources (compared to the no- wind case). With an infinite halo, the 
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probability density in the disk is given by 

dPt^L {rs\0} ^ f°° kJ Q (kr s )dk 

d 2 r s a i V c + Ky/V?/K 2 + 4fc 2 

K J_ Z" 00 xJ (a:) da; 

r« 7 r s /r wind + (r s / r winA ) 2 + a; 2 ' 

where the characteristic radius r W j n( j = 2K/V C has been defined. The expression in 
Eq. i|15|l is a function of r s /r w i n( j only. The deviation from a pure l/r s law, as well 
as deviations due to escape, radioactive decay and spallation (see next section), is shown 
in Fig. El The nj m -probability is given by 




Fig. 5. Deviation from the pure l/r s density profile, N(r s )/(l/4nKr s ), due to the var- 
ious effects studied here: escape from the z — ±L boundaries, spallations and Galactic 
wind. In this latter case, the choice r sca i c = 2r W j n d has been made to show the similar 
behavior at large r s . The case of a radioactive species has also been shown. It should 
be noticed, however, that in most interesting cases, the scale length i rac j is much smaller 
than the others, so that in this case the propagation is dominated by radioactive decay 
and spallations and Galactic wind can be safely discarded. 



12 Taillet & Maurin: Spatial Origin of Cosmic Rays in Diffusion Models 



Jb(a;) w j r Uro / , , a C 



4x 1 r w ; n d V ?^,ind r wind 



X' 



111 



Him /?"wind + \/ r iim/ r wmd + 



dx . (16) 



Some values are indicated in Tab. 0] and plotted in Fig. 

It is interesting to note that the effect of 2r win d is similar (though not rigo rously 
identi cal) to the effect of L (see Fig. [SJ. As a matter of fact, this was noticed by 



Jones 



1 1978) who studied the propagation properties in a dynamical halo and provided a very 
simple picture (along with a rigorous derivation) of the effect of the wind. Consider a 
particle initially located at a distance z from the disk. It takes a time tdiff ~ z 2 / K to 
diffuse back in the disk. In the meantime, convection sweeps the particle in a distance 
z w = V c <diff ~ V c z 2 j K. Both processes are in competition and the particle will not reach 
the disk if z w > z. This define an effective halo size L* w K/V c . This is our parameter 
''wind up to a factor 2. 



5.2. Pure spallation 

The Galactic disk contains interstellar gas mostly made of hydrogen. When Cosmic Rays 
cross the disk, they can interact with this gas. This interaction may result in a nuclear 
reaction (spallation) , leading to the destruction of the incoming particle and to the cre- 
ation of a different outgoing particle (secondary). We present two approaches to the 
problem of diffusion in presence of a spallative disk. When the halo is infinite in extent, 
the solution may be obtained by using the interpretation of diffusion in terms of random 
walks. This will be treated in App. C. In the general case, the Bessel developments can 
be used as before. Starting from Eq. (|A.4|I . the expression for the probability density is 
readily obtained. The limit L — > oo is noteworthy, as the resulting expression isolates the 
influence of spallations: 

dVT^ {r, |0} ^ f°° U (kr s ) ^ = 1 r°° xJ (x)dx 



d 2 r s J 2hT inc i + 2kK 4ir Kr s J r s /r spal + a; 

where the quantity r spa i = Kj (hTi ne \) has been defined. Would there be no spallation, the 
l/r s behavior would be recovered. The term 2hTi nc \ has the effect to kill the contributions 
of fc < fcgpai in the integral, with fc spa i = hT- me \/K. It leads to a decrease of the integral 
on scales r > r spa i = l/fc S pai- Some typical values, for K — (5KqR s (see Sec. I5.5fl with 
K = 0.03 kpc 2 Myr" 1 and <5 = 0.6 are given below at 1 GeV/nuc and 100 GcV/nuc. The 
heavy species are more sensitive to spallations, so that they come from a shorter distance. 
This could in principle affect the mean ato mic weight of Cosmic R ays if the composition 



of the sources is not homogeneous (see e.g. iMaurin et alJ l|2003a|) ). See Sec. 17. If) for the 
results with realistic propagation parameters. 
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p O Fe 


a (mb) 
r sp ai (kpc), 1 GeV/nuc 
7" sp ai (kpc), 100 GeV/nuc 


44 309 760 
10.2 1.45 0.59 
115 16.4 6.7 



Table 3. Some values of the inelastic cross section and the associated spallation scale 
length. 



For small values of r s /r spa i, the convergence of the previous integral is slow, and other 
forms obtained by integration by parts, as developed in the App. B.3, might be preferred. 
However, in this particular case, the identity 

/>oo 

xdxJ (x)/{x + a) = / ydye~ ay /(l + y 2 ) 3/2 



yields the more useful form 



d 2 r s iTvKr s J (1 + y 2 f/ 2 dy ' (18) 

This expression is in full agreement with Eq. (|G.3J| obtained with the random walk 
approach (see App. C). For large values of r s /r spa i, the convergence can be checked by 
comparing the results to the asymptotic development 
xJ Q (x) , 1 9 225 

— dx^ — J + — + ... 

x + a or or or 

Finally, the ru m - probability can be computed as before 

(< Him) = / ^ 3/2 1 ~ e ^}~ X ~} 

JO (1 + ST) ' L I r spal J 

Some values are indicated in Tab. 



5.3. Comparison and combination of the different effects 



To summarize, the effect of spallation and Galactic wind depends on the two parameters: 

(19) 



2K , / K(E) \ /lOkms" 

^wind = -Tj- ~ 5.87 kpc x ' ' * 



V c V0.03kpc 2 Myr"V V K 



K o 1 / Jif(S) \ / 100 mb 

r spal = — — w 3.17 kpc x - x rr 

"1 incl P \0.03 kpc z Myr"V \ o 

The ni m -probability is displayed in Fig. as a function of rn m . The effect of the Galactic 

wind is very similar to that of the top and bottom boundaries, whereas the effect of 

spallations is quite different. In the latter case, the cutoff in the density is a power law 

in r s and decreases much more slowly than the exponential cutoff due to the wind or to 

escape. As a result, the 99%-surfaces are much larger than the 90%-surfaces. This can 

also be seen in the first three lines of Tab. 0] 

When all the effects above are considered, Ea. l|A.4(l gives 



dV{r s \0} 1 f°° xJ (x)dx 
cx ' 



d 2 r s kr s 



Pspal + Pwind + V^wind + ^ COth { VPwind + ^ I PL } 
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Fig. 6. Integrated probability that a particle detected at the origin was emitted inside the 
ring of radius r s , in the three situations considered. The solid dark line is obtained when 
only the leakage through the z — ±L boundaries is considered, in which case the radii 
scale as r/r sca \ c with r sca i c = L. The dotted, respectively dashed, line is obtained when 
only the spallations, respectively only the convective wind, are considered. The solid grey 
line indicates the probability that the primary progenitor of a secondary detected in the 
Solar neighborhood was emitted from within a given distance. 

where p spal = r/r 

spah Pwind = r/r wind et p L = r/L. The smallest of these three num- 
bers indicates the dominant effect. Various rn m -probabilities are shown in Tab. 01 For 
a radioactive species, the spallations and the Galactic wind have a negligible effect on 
propagation as long as l ra( j (see Sec. 14. 2f) is smaller than L, r spa i and r W i n d- 

5.4. The number of disk-crossings in the general case 

Several properties (energy losses, amount of reacceleration, secondary-to-primary ratio) 
of the Cosmic Ray flux detected in the Solar neighborhood are determined by the number 
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L(kpc) 


Vc (km s" 1 ) 


<7i llc i (mb) 


Rso% (kpc) 


i? 90 % (kpc) 


#99% (kpc) 


Lio x 10 








6.3 xLio 


19 xLio 


34 xLk, 


oo 


10 x Vio 





4.9/Vio 


18.6/Vio 


41/Vio 


CO 





100 X CTioo 


3.7/aioo 


30.8/aioo 


318/ctioo 


5 








3.1 


9.5 


17 


5 


10 


50 


2.05 


6.8 


13.1 



Table 4. This table indicates the radius of several /-surfaces, for several values of L, 
V c and (Jinel- We have introduced L\o = L/IQ kpc, Vio = V c /10 km s _1 and cioo = 
cr/100 mb. 



of times a given Cosmic Ray has crossed the disk since it was created. The distribution 
of disk-crossings is computed in App. C in the case of an infinite diffusive volume and 
in the absence of Galactic wind. In the most general situation, the mean number of 
crossings (though not the entire distribution of crossing numbers) can be computed as 
follows. Each time a particle crosses the disk, it has a probability p = 2ft,(Ti nc inisM of 
being destroyed by a spallation. The number N(r) of surviving particles can thus be 
obtained from No(r), the number of particles diffusing without spallations, as 

N(r) =N (r) x — , 

so that the number of crossing is readily obtained from the densities with and without 
spallations as 

Notice that this expression applied to Eq. I|18H leads to Eq. (|C.2|I when L — * oo, V c = 0, 
and when p is small. As the surface density of the disk is 2/misN ~ 10~ 3 g cm -2 , the 
mean column density crossed by the particle (called grammage) is given by 



S(r s ) = n cmss (r s ) x 2hni SM ~ 20 g cm x 



\l Across (^"s) 



10 4 



The evolution of the grammage with the distance of the source is displayed in Fig. [7\ 
The effect of escape, spallations and Galactic wind is shown. 

As a cross-check, it can be noticed that in this approach, the mean grammage 



(S) S p atia i = / S(r.) 5 d 2 r s 



yields the right order of magnitude for the usual grammage derived from leaky box 
analysis (~ 9 g cm -2 ). Moreover, the knowledge of n cross (r s ) allows to estimate the 
magnitude of energy losses and reacceleration rates. 
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Fig. 7. Grammage crossed as a function of the origin, for some of the models discussed 
in the text and for a typical value of K = 0.03 kpc 2 Myr -1 . 

5.5. The energy dependence 

The diffusion coeffic ient actually depends on energy. A commonly used form (see 



Maurin et al 



1 2002b) for a discussion) is 



*' =a '"Sgv, 

where 1Z stands for the rigidity, Kq ~ 0.01 — 0.1 kpc 2 Myr -1 and 6 ~ 0.3 — 1. The previous 
results were given for K — 0.03 kpc 2 Myr -1 , typical for a proton with an energy of 1 GeV. 
This implies that the parameters r w ; n d , r spa i are larger at higher energy. They eventually 
become larger than L, so that at high energy escape dominates. At low energy, the relative 
importance of spallation and convection can be evaluated by comparing r w j n( j and r spa i. 
However, it must be noticed that even when r w j n( j is greater than r spa i, the Galactic wind 
may have a non negligible effect on the Cosmic Ray spatial origin because the cutoff due 
to r W ind is much sharper (see Fig. EJ. Moreover, the influence of the Galactic wind on 
the spectra is important because of the induced energy changes (adiabatic losses). 
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For the sake of dcfinitcness, we will consider from now on that the cosmic ray sources for 
stable primaries are located in the disk and that their radial distribution w(r s ) follows 
that of the pulsars and supernovae remnants, given by 



wsN(r s ) 



exp 



-/3 x 



R, 



(21) 



with Rq = 8.5 kpc, a = 2, /3 = 3.53 for 



Case fc B hattacharva (1998). This distribution is 



now closer to the distribution adopted by 



Strong fc Moskalenkc 



(a = 0.5 and j3 



1), a flatter distribution designed to reproduced radial 7-ray observations (see Fig. Ilfijl . 
This distribution can be inserted in Eq. 0] which is then used to compute the /-surfaces. 
These surface are displayed in Fig. [S] for three cases (L = 2 kpc, L = 5 kpc and L = 
10 kpc). For large halos, the source distribution acts as a cutoff and greatly limits the 
contributions from peripheric Galactic sources. 

The results are not much affected by taking an a ngular depen dence into account. 
Considering for example the spiral arms modelling of Vallcc (2002), Fig. [S] shows that 
the extension of the /-surfaces is almost not affected by these small scale structures. In 
the rest of this paper, the purely radial distribution <|21[) is assumed. 
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Fig. 8. 99%-surfaces for R = 20 kpc and three cases, L = 2 kpc, L = 5 kpc and 
L = 10 kpc. 
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7. Application to the propagation parameters deduced from the observed 
B/C ratio 

The previous sections present a complete description of the origin of Cosmic Rays in a 
stationary diffusion model (energy losses and gains are discarded). To each process by 
which a Cosmic Ray may disappear before it reaches the Solar neighborhood is associated 
a parameter: L (escape through the top and bottom boundaries), r w ; n d (convection), r spa i 
(destructive spallation). The relative importance of these parameters may be measured 
by the two quantities Xwind = -k/fwind an d Xspai = L/r spa x- One can distinguish three 
regimes which determine the diffusion properties of the system: i) the escape through 
the boundaries dominates for Xwind <C 1 and Xspai ^ 1; ii) convection dominates for 
Xwind > 1 and Xwind > Xspai; hi) spallations dominate for x spa i > 1 and Xwind < Xspai- 
We now use the s ets of diffusion parameters consistent with the B/C data given in 
1 2002af) (hereafter MTD02) to evaluate realistic values for these quantities. 



Maurin et al 



7.1. Evolution of Xwind and Xspai with S 

In MTD02, we provide for each configuration a (source spectral index), 5 (diffusion 
spectral index) and L (diffusive halo size) the corresponding Ko, V c and V a (Alfvcnic 
wind responsible for reacceleration) that fit best the ratio B/C. In this study, V a is not 
very important since it only changes the energy of the particles: a Cosmic Ray emitted at 
1 GeV/nuc and gaining a few hundreds of MeV/nuc during propagation will be detected 
at a slightly larger energy, for which the results given here will not be very different. This 
becomes even more true beyond a few GeV. Reacceleration will be ignored throughout 
this study, as well as energy losses, for the same reason. Moreover, the values of Ko, V c 
and V a do not depend much on a (see Fig. 9 of MTD02), so that Xwind and Xspai depend 
mainly on S and L. They depend on rigidity, through K(E), as can be seen in Fig. [9] 
where Xwind and Xspai are displayed as a function of 5 for several species, several values 
of L and several rigidities. 

The left panel displays Xwind {$>L) f° r three rigidities: 1 GV, 10 GV and 100 GV. 
Up to several tens of GV, convection is in competition with escape; afterwards escape 
dominates. The noticeable fact is that models corresponding to S < 0.45 are escape- 
dominated, whereas convection dominates only for large S at low energy. It appears that 
all other parameters being constant, Xwind is fairly independent of L. (indicating a similar 
relative importance of convection and escape for the models reproducing the B/C ratio, 
see MTD02). However, the spatial origin does depend on L and 7- w ind and not only on 
their ratio. 

The right panel of Fig. [5] plots Xspai(<5, L) for 1Z = 100 GV and 1 GV for various 
nuclei. Protons are the most abundant species in Cosmic Rays. Boron and CNO family 
are important because they allow to constrain the value of the propagation parameters, 
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Fig. 9. Left panel: Xwind(5, TZ) as a function of the diffusion spectral index S for different 
rigidities TZ; from top to bottom, TZ = 100 GV, TZ = 10 GV and TZ = 1 GV. The 
parameter Xwind, as well as Xspai, is not very sensitive to the halo size L. Right panel: 
Xspai(<5, TV) as a function of S for TZ = 100 GV (upper curves) and TZ = 1 GV (lower 
curves) for four species: p [a ~ 40 mb), d (a ~ 100 mb), B-CNO (ct ~ 250 mb) and Fe 
(cr ~ 700 mb). For the latter species we plotted the same three L values as in left panel. 
The behavior for other species is similar so that we only plotted the case L = 6 kpc. 



e.g. through the B/C ratio. Last, the Fe group provides another test of the secondary 
production via the sub-Fe/Fe ratio. The evolution of x spa i for these species is conform 
to what is very well known from earlier leaky box inspired studies: for heavier nuclei, 
spallation dominates over escape and for this reason, the induced secondary production 
is particularly sensitive to the low end of the grammage distribution. 

To summarize, the left panel shows the evolution from convection-domination to 
escape-domination as a function of TZ and 5, the effect of the wind being negligible 
above ~ 100 GcV whatever S (xwind ^ 10). The right panel gives the evolution from 
spallation-domination to escape-domination as a function of TZ, S and the species under 
consideration. The effect of spallation is more important for heavy than for light nuclei, 
but this difference is too small to prod uce an evolution of t he average logarithmic mass 



for high energy (~ TeV) Cosmic Rays l|Maurin et al 



2003a). 



7.2. Spatial origin in realistic diffusion models at 1 GeV/nuc 

From the previous discussion, it appears that spallations and Galactic wind play a role 
at low energy. The results will be shown for the particular value 1 GeV/nuc which is 
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interesting for various astro-physical problems. First, once modulated, it corresponds to 
about the very lowest energy at which experimental set-ups have measured Galactic 
Cosmic Rays. Second, the low energy domain is the most favorable window to observe 
p (resp. d) from exotic sources (see companion paper iMaurin fc Ta illet (2002)), as the 
background corresponding to secondaries p (resp. d) is reduced. Last, these energies 
correspond to that of t he enduring probl e m of the diffuse GeV 7-ray radial distribution . 



This was first quoted by 



Stecker fc .Tone; 



■Tone; 



1 1977) and further investigated by 
taking into account the effect of a Galactic wind. 

From the sets of diffusion parameters that fit the B/C ratio, the values of the pa- 
rameters r w j n( j and r spa i are computed (see Tab. |SJ for the four nuclei shown in Fig. [5] 
and for three values of 5 = [0.35, 0.6, 0.85]. From these values, the 50-90-99%-surfaces are 





P, P 

S = 0.35/0.6/0.85 


B-CNO 


Sub-Fe, Fe 


L — 10 kpc r wind (kpc) 
r spa i (kpc) 
\J (r 2 )(kpc) 


co/5.17/1.6 
33.5/10.2/4.0 
6.43/4.67/2.50 


oo/3.41/0.89 
4.21/1.07/0.35 
5.23/2.92/1.11 


co/3.26/0.83 
1.46/0.37/0.12 
4.00/2.03/0.57 


L = 6 kpc r wind (kpc) 
r spa i (kpc) 
v/ (r 2 ) (kpc) 


oo/3.64/1.15 
24.7/7.4/2.9 
4.93/3.5/1.88 


oo/2.40/0.64 
3.10/0.80/0.26 
3.96/2.18/0.82 


oo/2.30/0.60 
1.08/0.27/0.09 
2.99/1.63/0.54 


L — 2 kpc r-wind (kpc) 
r sp ai (kpc) 
yj (r 2 )(kpc) 


co/1.40/0.46 
9.7/2.9/1.2 
2.07/1.44/0.78 


oo/0.92/0.26 
1.21/0.31/0.10 
1.63/0.87/0.33 


co/0.88/0.24 
0.42/0.10/0.03 
1.21/0.57/0.19 



Table 5. Values of r w i n d and r spa i for the sets of parameters that, for a given 5, give the 
best fit to the observed B/C ratio. The mean square value yj (r 2 ) of the distance to the 
sources is also shown. 



derived and displayed in Fig.^J for protons and Fe nuclei. The effect of 5 fFig. H2|) . of L 
(Fig. I13f) and of the species fFig. Ill[). are considered separately. As regards the first two 
effects, the /-surfaces are smaller: (i) for greater values of <5, mainly because the effect 
of the wind is then greater, and (ii) for small values of L, as in this case escape is more 
important. 

As regards the last effect, it can first be seen from Fig.ll3lthat the heavier species come 
from a shorter distance (because the spallations are more important). The secondary 
species can be treated simply by using a source function obtained by multuplying the 
primary density by the gas density. It would be straightforward to apply the previous 
techniques to a realistic gas distribution (taking into account, in addition to the fairly 
flat HI distribution, that of molecular H2 and ioniz ed HII which are more strongly peaked 



in the inner parts, see e.g. 



Strong fc Moskalenkol l|1998fl for a summary and references) 
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and to infer the contours inside which the secondaries are created. The corresponding 
/-surfaces are not shown here, as they would be quite similar to those of the primaries 
(see left panel) . What we do display in the right panel are the /-surfaces of the primaries 
that lead to given secondaries, as these progenitors determine the secondary-to-primary 
ratios (see Sec. 14.10 . 




Fig. 10. (50-90-99)%-surfaces (protons and Fe nuclei are considered), in a typical diffu- 
sion model with L = 6 kpc and 5 — 0.6. 




Fig. 11. 99%-surfaces for several species. The left panel corresponds to primary species 
(protons, CNO and Fe) while the right panel corresponds to the progenitors of secondary 
species (B and sub-Fe), for L = 6 kpc and S = 0.6. 

7.3. Effective number of sources 

From the previous results, it appears that only a fraction of the sources present in the disk 
actually contribute to the flux in the Solar neighborhood. In this paragraph we present 
the fraction f s of the sources which are located inside given /-surfaces. This fraction is 
presented in Fig. El for the particular model L — 6 kpc, S = 0.6, for protons and Fe 
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Fig. 12. 99%-surfaces for several 5, in the case L — 6 kpc. The left panel corresponds to 
protons while the right panel corresponds to Fe nuclei. 




Fig. 13. 99%-surfaces for several L, in the case S = 0.6. The left panel corresponds to 
protons while the right panel corresponds Fe nuclei. 

nuclei. It can be read that for example, it takes 7.6% (resp. 1.5%) of the Galactic sources 
to make 90% of the protons (resp. Fe nuclei) reaching the Solar neighborhood. 

Figure 1141 also shows that a very small fraction of sources may contribute to a non 
negligible fraction of the fluxes. For example, the sphere of radius r ~ 100 pc centered 
on the Solar neighborhood contains only 2.5 10~ 5 of the sources but for L = 6 kpc and 
5 = 0.6, it is responsible for about 5% of the proton flux and 18% of the Fe flux. The 
mean age of the Cosmic Rays is given by (t) ~ (r 2 ) /2K ~ 7-400 Myr (see Tab. [J. For a 
supernova rate of three per century, the total number of sources responsible for the flux 
is ~ 2 10 5 — 10 7 . This tells us that in models with the largest 5, 18 % of the Fe flux can 
be due to only 5 sources. The approximation of stationarity and of continuous source 
distribution is likely to break down with such a small number of sources. Conversely, for 
small values of 6 (preferred by many authors), this approximation is probably better. 
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Fig. 14. Fraction f s (in %) of the Galactic sources contributing to the fraction / of 
the Cosmic Ray flux at the Solar position, for protons and Fe nuclei, for the particular 
diffusion model L = 6 kpc, S = 0.6. 





P, P 

8 = 0.35/0.6/0.85 


Fe 


L = 10 kpc 90% 
99% 


22.9/13.3/3.9 
56.7/40.7/16.4 


9.6/2.2/0.24 
35/14.5/1.77 


L — 6 kpc 90% 
99% 


14.3/7.6/2.3 
40.9/27/9.6 


5.63/1.5/0.17 
22.2/9.2/1.48 


L — 2 kpc 90% 
99% 


2.9/1.6/0.83 
9.3/5.9/3.8 


1.2/0.24/0.025 
4.7/1.77/0.27 



Table 6. Fraction / s (in %) of the Galactic sources contributing to a given fraction (90% 
and 99%) of the Cosmic Ray flux at the Solar position, for protons and Fe nuclei, for the 
diffusion models studied before. 



7.4. Radioactive species, e + and e and the local bubble 



Donato et al 



1 2002) emphasized that the existence of a local underdensity (n < 
0.005 cm -3 ) around the Solar neighborhood greatly affects the interpretation of the 
flux of radioactive species at low energy (we refer the interested reader to this paper for 
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p, p Sub-Fe, Fe 

O — U.OUj/ u.u/u.oo 


L = 10 kpc 


196 / 392 / 332 76 / 74 / 17 


L = 2 kpc 


72/ 130 / 108 24/ 22 /7 



Table 7. Mean Cosmic Ray age (t) ~ (r 2 )/2K in Myr for the models studied in this 
paper. 



a deeper discussion and references on the local interstellar medium). The most impor- 
tant effect of this hole is that it exponentially decreases the flux by a factor exp(— a/i ra( j) 
(a < 65 — 250 pc is the radius of the local underdense bubble and Z ra( j is given by 
Eq. (|13[0 . This can be easily understood as there is almost no gas in this region, hence 
no spallations, leading to no secondary production. The local bubble is obviously not 
spherical, but this approximation is sufficient at this level. This attenuation factor is 
straightforwardly recovered starting from the probability density as given in Sec. 14.21 if 
correctly normalized to unity. To this end, the sources (here spallation of primaries on 
the interstellar medium) are considered to be uniformly distributed in the disk, except 
in the empty region r < a. The probability density is zero in the hole whereas outside, it 
is given by 

^ - exp ( JL) dV^ = 6Xp( : (r * ~ a)/l - d) rf»r. . (22) 

\trad / £7tTs - ^rad 

The quantity V^^{r s < ry im \0) is obtained directly from the no hole case (see Sec. 14.2(1 
by replacing ri; m by rg m = rn m + a. It means that the sources that contribute to the 
fraction / = (50 — 90 — 99)% of the flux of the radioactive species are located between 
a and rg m = (0.7 — 2.3 — 4.6) x ^ ra d + a. Hence, the hole only plays a marginal role for 
the origin of a radioactive species (unless ^ ra d < a )? whereas the result for the flux is 
dramatically different. 

We saw in a previous section that the high energy e + and e~ behave like unstable 
species. Their typical length ri oss can be compared to i ra( j 



Wkpc) _ Q 12 / n ^ I K (n/i cvy 



y/ro/l Myr ' V 1 GV ]j 0.03 kpc 2 Myr" 1 ' ^ 

n_(kp.)-^ x / A '" W1 GV)< 



E y 0.03 kpc 2 Myr^ 1 
The dependence in the propagation model is similar for both expressions and is contained 
in the last term. There is a big difference, though, as the typical distances travelled by 
radioactive nuclei scale as VlZ, whereas they scale as 1/VE for electrons and positrons. 

Realistic values for ? ra( j and n oss are presented in Fig. 1151 At high energy, the Lorentz 
factor enhances the lifetime of radioactive nuclei, making their origin less local, whereas 
the energy losses are increased for electrons and positrons, making their origin more 
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local (99-90-50% of 100 GeV e + e _ come from sources located in a thin disk with radius 
r fim e ~ ~ L1 ~ °- 55 ~ °- 38 k P c )- For 100 GeV e+ and a11 models fitting B/C have the 
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E k/nuc (GeV/nuc) E (GeV) 



Fig. 15. Realistic values of / r ad/\/iT)7lMyr and ri oss for two extreme halo sizes L and 
diffusion slope 8. As all results in this section, propagation parameters fit B/C and are 
taken from MTD02. 



same value for Kq1Z s , because at 100 GeV/nuc, spallations and convection are negligible 



I Ma urin et al 



2002a). As a consequence, for GeV energies, rioss increases more rapidly 



for small 6 than it do es for larger 8. To stud y the effect of local contributions to the 



spectra of e + and e , 



Aharonian et al 



1 1995|) used the value 6 = 0.6 and compared to 
other works with 5 = 0. As these authors noticed, the modelling in the whole energy 
spectrum is 6 dependent, but our study gives the range compatible with B/C studies. 

Finally, radioactive nuclei are a very important tool for Cosmic Ray physics. They 
come from a few hundreds of parsec, and their fluxes are very sensitive to the presence of a 
local underdense bubble, through the attenuation factor k = exp(— a/Z ra d)- For example, 
for a typical bubble of size a = 100 pc and an energy 800 MeV/nuc (interstellar energy), 
K o.35 ~ exp(— 0.33-y/l Myr/ro) if S = 0.35, whereas kq.85 ~ exp(— y^I Myr/r ). With 



- 10Bc = 2.17 Myr, r 26A1 = 1.31 Myr and t^ vi = 0.443 Myr, it leads to Kg 



0.80-0.75-0.61 and Ko. c 8 ' 5 A1, C1 « 0.51-0.42-0.22. For 14 C, the attenuation is k 14C < 1 
around 1 GeV/nuc, so that this species is heavily suppressed. However, it should be 
present around 10-100 GeV/nuc (as k 14C ~ 1 at these energies), with the advantage that 
solar modulation is less important at these energies. 

The flux of radioactive species directly characterizes the local diffusion coefficient Kq 
if the local environment is specified. This would in turn allow to break the degeneracy 
in propagation parameters that one can not avoid at present. Last, even though the 
surviving fraction of a radioactive does depend on the halo size L, we emphasize that it 
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is a very indirect way to derive the propagation parameters. In the forthcoming years, 
new measurements of radioactive species that do not depend on L (e.g. by PAMELA and 
AMS) should provide a promising path to update our vision of cosmic ray propagation. 



8. Summary, conclusions and perspectives 



The question of the source distribution is very present in Cosmic Ray physics. With the 
occurrence of the old problem of sh ort pathlengths dist r ibutio n in leaky box models (see 



for example 



Webber et al 



(1998)), 



Lezniak fc Webber 



19791) st udied a diffusion model 



Webberl lll993aIbD propagated 8 



with no-near source in the Solar neighborhood. Later, 
like sources with dif f usion generated by a Monte Carlo random walk for the same purpose. 
Brunetti fe Coding Ij2000(l follow this line but they introduce random walks in a more 
realistic environment, i.e. c i rcular , elliptical and spiral magnetic field configurations. In 
a more formal context. iLed (J1979) used a statistical treatment of means and fluctuations 
(see references therein) to characterize amongst others the possibility that nearby recent 
sources may dominate the flux of primaries. Finally, it is known that the present Cosmic 
Ray models are not able to reproduce accurately for example proton-induced 7-rays 
measurements. To illustrate this point, we plot in Fig. llGl thc radial distribution of protons 
obtained with the same diffusion parameters as used above. None of the models shown 
match the data. One is left with two alternatives: either modi fying the source distribution 
(for example, the distribution of lStrong fc Moskalen ko ( 1998) yields a better agreement), 
or giving up the assumptio n that the diffusion parameters apply to the Galaxy as a whole 



( Brci tschwerdt et al 



2002). It is thus of importance to understand to what extent the 



Cosmic Rays detected on Earth are representative of the distribution of the sources in 
the whole Galaxy. 

We provide an answer to this question under the two important hypotheses that the 
source distribution is continuous and that we have reached a stationary regime: most of 
the Cosmic Rays that reach the Solar neighborhood were emitted from sources located 
in a rather small region of the Galactic disk, centered on our position. The quantitative 
meaning of "rather small" depends on the species as well as on the values of the diffusion 
parameters. For the generic values <5 = 0.6 and L = 6 kpc chosen among the preferred 
values fitting B/C (see Sec. 0), half of the protons come from sources nearer than 2 
kpc, while half of the Fe nuclei come from sources nearer than 500 pc. Another way to 
present this result is to say that the fraction of the whole Galactic source distribution 
that actually contributes to the Solar neighborhood Cosmic Ray flux can be rather small. 
For the generic model just considered, 8 % (resp. 1.5 %) of the sources are required to 
account for 90 % of the proton (resp. Fe) flux. These fractions are smaller for higher <5 
and smaller L. To summarize, the observed Cosmic Ray primary composition may be 
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radius (kpc) 



Fig. 16. Radial distribution of the proton flux for the models di s cusse d in this study, 



compared to the source radial distribution of 
For each of the values L — 2 kpc and L 



Case k Bhattacharval jl99jj) given Eq. J2TJ- 
10 kpc, the three values 8 = 0.35, 



0.6 and 0.85 are presented, the flatter distribution corr esponding to the lower 5. Also 
shown is the gamma-ray emissivity per gas atom (COS-B iBloemerj J 1989 ), which is pro- 



Bloemenl |1989)) and 



portional to the proton flux, as given by COS-B (open circles, 
eg ret (triangles, IStomg fc Mattoxl Ijl996l^ . along with the proton flux obtained with 



the 



Strong fc Moskalenkc 



distribution (see Sec.EJ). 



dominated by sources within a few kpc, so that a 



particular care shou ld be taken to 



2003bh. 



model these source, spatially as well as temporally IjMaurin et al. 

Independently of all the results, this study could be used as a check for more so- 
phisticated Monte Carlo simulations that will certainly be developed in the future to 
explore inhomogeneous situations. Several other consequences deserve attention. First, 
the results may point towards the necessity to go beyond the approximations of both 
continuity and stationarity. In particular, it could be that only a dynamical model, with 
an accurate spatio-temporal description of the nearby sources, provides a correct frame- 
work to understand the propagation of Galactic Cosmic Rays. The contribution from 
nearby sources would be very different in the low energy (~ GeV/nuc) or in the high 



energy regime (~ PeV) compared to the stationary background. Second, as discussed in 
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Sec. 14.11 the diffusion parameters derived from the observed B/C ratio have only a local 
validity, and one should be careful before applying them to the whole Galaxy, since the 
Cosmic Rays are blind to most of it. 
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Appendix A: General solutions of the diffusion equation in cylindrical 
geometry 

A.l. General solution for our diffusion-convection model 

For a primary species, t he differential density (in energy) N(r, z) is a solution of the 



equation (see for example iMaurin et al.l l)2002b|) and references therein) 



£ di &N(r, z) + r rad AT(r, z) + 2h5(z) n 1SM .cr incl .v N{r, z) = 2h8(z)w(r) , (A.l) 

with £,jiff = — K + rir( r 3r)) ■ ^ ne various terms in Eq. (|A. 1|) correspond 

respectively to (i) a differential operator £diff describing convection V c out from the 
Galactic plane and isotropic diffusion K throughout the confinement volume; (ii) ra- 
dioactive decay of the unstable nucleus in the whole Galaxy; (iii) destruction in flight 
n\sM-^mei-v when crossing the gaseous thin disk of constant height 2h and constant 
density nisM! (iv) a source term in the thin disk. The solution is determined by the 
boundary conditions of the problem where Cosmic Rays freely escape, i.e. we demand 
N(r — R,z) — N(r, \z\ = L) = (L is the half-height of the diffusive halo, R the radial 
extension of the Galaxy) . The solution in the disk (z — 0) is given by 



%0) = Ef Jb(C^) (A.2) 



Q 

with 



2hT incl + V C + KSi coth and Sf = ^ + ^ + 4^ . (A.3) 

and Qi is the Bessel transform of the source distribution (which may depend on the 
density of another species, in particular for secondary species). 

For a primary point source, w{r) = 8(r) / (2irr) and we find in the disk {z = 0) 



^V,0)=E{^ WI :}xJ (c4). (A.4) 
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The generi c solution for second aries can be straightforwardly derived from that of 

(20oJ), 



primaries (e.g. 



Maurin et al 



N c y\r,0) = 2ftr pri „ c x]T 



1 



i l7rJ 1 2 (c i )i^ 2 A p^l^, ^ 



( Q i ) • 



(A.5) 



The distinction between Af rim and A s - cc is neces- 



We USe Tprim^scc = niSM^^prinnscc- 

sary since both species have different destruction rates and rigidities. 



Appendix B: Numerical evaluation of the point source solution in Bessel 
basis 

In practice, the infinite sums above are truncated to some order ntronc, chosen as a 
compromise between accuracy (good convergence of the series) and computer time. In the 
case of a point source S(r), the profiles are singular near the source and the convergence 
of the series appears to be very slow. A few methods are presented to speed up this 
convergence. 



B.l. Softening of the source term 

First, the source term may be spread out on a radius a, by replacing the ^-function by 
w(r) = 9(a — r) /(ira 2 ) for which an extra 2Ji(Qa/R)/(Qa/R) term appears in the Bessel 
transform. With a judicious choice of the parameter a, the solution is very close to the 
original for r»a, but convergence is much faster due to the extra 1/Q factor. 



B.2. Sum representation: comparison to a known function 

Part of the difficulty to evaluate numerically the Bessel expansions comes from the fact 
that the resulting functions are singular at the source position. If we know a reference 
function / ref (r) which exhibits the same singularity and for which the Bessel coefficients 
1^1 are known, it is then judicious to write the density (p=r/R) as 

+ / rof (p)- (B.l) 

where the singularity is entirely contained in the / term, so that the Bessel expansion 
has been regularized. The choice of / may be guided by the behaviors observed in Sec. [31 
For L ~ R, the solutions should be quite similar to the solution in spherical geometry, 
given by / rcf (r) = (1 — p)/(A-KpKR). For very small halo size [L < 1 kpc), the effects 
of the top and bottom boundaries are dominant and a modified function / (r) = (1 — 
p)j (ApirKR) x exp(~pR/L) is more adapted (see Eq. ©). This method yields a very 
good and rapid convergence as long as sources are located in the thin disk z = 0. 



Nf(r,0) 



"tronc 

E 

. i=l 
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B.3. Integral representation for infinite radius disk 

When the disk has an infinite radius, the Bessel sum can be replaced by an integral, 
and the end result is obtained from the Bessel sum by the substitution Q/R — * k and 
w ttQ/2 — > knR/2, so that in the general case - see Eq. <|A.4() and Eqs. (JA.2I) , 

nL , , f-V c z\ r kJ o(kr) sinh{5(fc)(L-z)/2} 

AT (r,^)=exp^— jy 2iM(fc) sinh {5 (fc)L/2} ' (R2) 



with 



A(fc) =2^r inc i + y c + X5(fc)coth^^^^ and S 2 (fc) = |% + + 4fc 2 . (B.3) 

The integrals of the form 

/>oo 

/[/] = / J (x)f{x)dx (B.4) 
Jo 

where f{x) is a function such that /(oo) = 1 are quite tricky to compute numerically, 
due to the very slow decrease of the oscillations in the integrand. Two remarks are of 
great help. First, as Jo(x)dx = 1, we have 

/[/] = 1 - / J (x) (1 - f(x))dx . (B.5) 



The convergence is faster, as 1 — / — > when x — > oo. Second, using the identity 
(xJi)' = xJq and integrating by parts, one has 



I[f\ = J~Mx) (^-f(x)^dx 



(B.6) 



This expression is meaningful only if xf(x) is a bounded function near the origin x = 0. 
Using the identity Jq = — J\ and integrating by parts again, 



J (x) [ f"{x) ~lM + lW) dx . (b.7) 



'[/]= , , 

\ x 

The latter expression is meaningful only if x 2 f(x) is a bounded function near the origin 
x = 0. These expressions provide several efficient alternatives to evaluate /[/]. 

Appendix C: Alternative description of spallations: random walk approach 

A Cosmic Ray crossing the Galactic disk has a probability p to disappear in a nuclear 
reaction with interstellar matter. This probability is related to the reaction cross section 
cr by 

a 



p = nicrni SM h = 6 x 10 Ki , - 

100 mb 

where K\ ~ 1 contains the dependence on the incidence angle of the particle with the 
Galactic plane. The propagation in the z axis is a one-dimensional random walk, so that 
for a Cosmic Ray emitted in the disk and reaching again the disk after a number t* of 
elementary random steps, the probability distribution of disk-crossing numbers n is given 
by (Papoulis 2002) 



dVMt, z(t) = 0) = exp dn . 

K 2 t* \ K 2 t* , 
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In this expression, t* is the number of steps of the walk z — 2j i=1 Zi and K2 ~ 1 depends 
on its statistical properties (for instance, K2 = 2 for elementary steps 2j = ±A and 
ft 2 ~ 1-43 for Zi uniformly distributed in the interval [—A, A]). The diffusion coefficient is 
defined as 

<z 2 ) K 3 X 2 t* K3 X 2 1 2 
A = = = = - K3TV , 

2t 2t 2r 2 

where k 3 = (z 2 ) /A 2 is the variance of the elementary random step (in units of A), so that 

the physical time is related to t* by t — t* x r = t* x 2K/k 3 v 2 , where A is the mean free 

path and v the velocity. We thus finally have, 

, , , , , , AKn ( 2Kn 2 \ , 

dV(n\t, z(t) = 0) = 5- exp 5- 

\ K2K3VHJ 

We are now able to compute the probability distribution of disk crossings for Cosmic 
Rays emitted from a distance r in the disk as 

dV{n\r,z(t)=0) = f 00 = y v , t < z{t) = Q) dt / c>1 n 

an J 

where the probability that a CR reaching distance r in the disk was emitted at time t is 
P(t|r,*(t)=0)«^^e«p v 

The above integral lC.il can be performed, yielding the final result 
dP(n\r,z(t)=0) r 2 n ' .2^2 \ -3/2 



> il+-V > ( C - 2 ) 

an 

with r 2 = 8K 2 / K2K3V 2 . The average number of disk crossings is readily obtained: 
dP(n\r,z(t) = 0) r 



n- 



■ dn = — 



dn ro 

and the associated variance is tends to infinity. We can also compute the integrated 

probability, that more that no crossings have occurred, as 

2 2\ -3/2 



T(n > n \r,z(t) = 0) = I 1 + 

A particle having crossed n times the disk has the probability p n = (1 — p) n ~ exp(— np) 
of surviving, so that the survival probability at distance r is given by 
[°° dP(n\r,z(t)=0) f°° xdx xrv/rn 

SUrv( )_ io ^ dn -J (i + ,2)3/2 e 

This can be written as 

V (r) - f°° p -otxr/r spa ,i 

Jo (1 + ar) ' 

with r spa i defined in Eq. (|19(l and a = Kiy K2K3/2. The density of Cosmic Rays in the 

disk is then given by 

1 1 Z" 00 t d-r 

N ( r ) = ~a — 17~ x = -— / _ /2 e-""/'-— . (C.3) 

AnKr AnKr J Q (1 + a; 2 ) 3 ' 2 

The quantity a seems to be related to the detailed statistical properties of the random 

walk under consideration, through K\, K2 and K3. However, direct comparison with the 

alternative expression l|18(l obtained above indicates that these expressions are indeed 

equivalent, with a = 1. 
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